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Using the Bardeen-Cooper-Schriefer (BCS) ansatz, we investigate the deviations in the behavior 
of ground state thermodynamic properties of an interacting homogenous Fermi gas from the zero 
range interaction approximation. The considered finite range interatomic potential allows us to vary 
the strength of the interaction to study the whole crossover from the BCS limit to the Bose-Einstein 
condensation (BEC) limit with no approximations. In particular, we concentrate in analyzing the 
thermodynamic Contact variable across the weak and strong coupling limits. Our analysis validates 
the correctness of the zero range approximation (contact approximation) in the strongly interacting 
regime. We predict quantum phase transition-like behaviour for long range interatomic potentials. 
These findings follow from an energy dependent gap due to the finite range of the interatomic 
potential. 



The experiments performed on dilute ultracold two- 
component Fermi gases represent perhaps the cleanest 
scenario that allows to access the crossover from Bardeen- 
Cooper-Schrieffer (BCS) superfluidity of Cooper pairs 
to Bose-Einstein condensation (BEC) of tightly bound 
fermion pairs [iKH]. These systems can be appropri- 
ately modelled by a Hamiltonian describing interact- 
ing fermions in a homogeneous environment. Inhomo- 
geneities created by a confining potential may be in- 
cluded via the local density approximation (LDA) [7]. 
Although the study of the crossover has been extensively 
addressed in the literature at zero [8H^ and non-zero 
temperatures |26ll27j . still, important issues remain to be 
discussed. One of them refers to the recently identified 
contact variable [2H] ? which besides from being amenable 
to direct measurement, yields information on a host of 
different properties, among them, density-density corre- 
lations 29 , the number of pairs in the closed channel of 
a Feshbach resonance [301 , and on the forms of ther- 
modynamic variables at unitarity |28| . Additionally, as 
we argue below, the calculation of this quantity is very 
sensitive to different approximations and schemes. In 
this Letter, we readdress the crossover within the BCS 
mean-field theory (MF) but without making the usual 
contact approximation in the interatomic interaction po- 
tential, C/ (r) « (5 (r) ^'^^ " where a is the s-wave scat- 
tering length. Namely, we calculate the thermodynamics 
at zero temperature for a realistic finite range potential. 
Although there have been previous analysis with finite- 
range interatomic potentials [311 [33], the issues addressed 
here have not been discussed. 

The system is a balanced mixture of two monoatomic 
fermionic species with the same atomic mass m, and 
whose interatomic potentials is U{\r — r'\). The Hamil- 
tonian in momentum space is given by 
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where V is the volume that confines the sample, e^, — 
h?k'^ I2m is the energy of a particle with momentum ±/ifc. 
flfc and bk are the annihilation operators of the two dif- 
ferent species and obey anticommutation rules. Uq is 
the Fourier transform of the interatomic potential U {r) . 
Here, for simplicity we have assumed that the interatomic 
potential is the same for any pair of fermions. 

The BCS ansatz [3iH35] for the ground state of the gas 
assumes that this state is a superposition of all possible 

Cooper pairs in the fiuid, |*bcs) = Ilfc (""fe +^fc^l^-fe) |0) 
being Uk and variational parameters with the prop- 
erty u^. A- vf, = 1 ensuring the normalization of the state. 
A straightforward evaluation in the grand canonical en- 
semble of the expectation value of the grand potential 
VL = {^^c?,\H - MEfc(4afe + blKW^cs), where the 
common chemical potential /i for A and B species en- 
sures that the mixture is balanced, yields 

17 = 2^(6, - er )vl + ^ E Uk-k'FuFu-: (2) 
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where Fk — UkVt- The term stands for the Hartree 
and Fock energy contributions, given by, 

^r-|E4-2^E^---^ (3) 

A straightforward variational scheme requiring that 
d^l/dvk = for all k, yields the well known equation 
for the gap, 

Ak^--YUk-k'—, (4) 

k' 

where the quasiparticle excitation energy is — [(efe — 
fi + f-k^Y + ^fe]^^^- supplements the above equa- 

tions with the total number of particles N = 2^j,t;^. 
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The above equations are valid for any short range in- 
teratomic potential. The contact approximation |36| is 
implemented by setting all the momenta dependence as 



and as it is well known, such a proce- 



dure inevitably leads to convergence difficulties at large 
values of the momentum q, with the necessity of either, 
schemes of renormalization nature or introducing physi- 
cally justified cut-ofl" parameters (3^1 EZI- This approx- 
imation disregards the significant role of the range a 
of the interatomic potential in determining the physics 
of the problem. In order to avoid this approximation 
we use a physically acceptable interatomic potential, an 
exponential-type of potential. 
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where T^o is the depth of the potential, and whose Fourier 
transform is 



(l + fcV2)' 
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This potential has essentially two main virtues. First, as 
a realistic potential it depends on at least two parame- 
ters, Vq and a. Second, but more important for our pur- 
poses, as derived by Rarita [38], the two-body scattering 
length is exactly known, 
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where x = a^/IVQmjh^ Jo{^) a-nd No{x) are zero order 
Bessel functions of first and second kind respectively with 
7, the Euler-Mascheroni constant. Therefore, for fixed 
range a, the Grand Potential O depends on the scatter- 
ing length via Vq. As Vq varies, a passes through an 
infinite series of scattering resonances where a diverges 
changing abruptly its sign; the resonances occur as ac- 
tual bound states emerge [38 . As it has also been widely 
discussed, those potential resonances are good approxi- 
mations to experiments since essentially all crossover ex- 
periments are in the so-called broad Feshbach resonance 
limit where the width of the resonance is much larger 
than the Fermi energy [29l |33l |39] . In this study we con- 
sider the BCS-BEC crossover through the first resonance. 

To analyze the thermodynamics in the whole BCS- 
BEC crossover, namely, as the inverse scattering length 
varies from — oo to c», we numerically solve [40j the set 
of equations Q-Q without any further assumption or 
approximation. Although not an independent thermo- 
dynamic variable, we also discuss the behavior of the en- 
ergy dependent gap A^. Unlike the standard BCS-BEC 
crossover mean-field theory by Leggett, the use of a finite 
range potential causes this energy dependence. We de- 
termine the chemical potential /i, the pressure p = — fl/V 
and the contact variable C given by [28l |30l |4T] 



C = - 
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FIG. 1. (color online) (a) Dimensionless gap A/Amax as a 
function of interatomic interaction range a and as a function 
of the scaled energy ecr^, with e = ek/ep- (b) Quasiparticle 
excitation energy, E = jl — Ej., for different values of the 
potential range a. The thick line corresponds to k. These 
calculations correspond to = 0. 



All variables can be adimensionalized in terms of the 

h k 

Fermi energy and momentum, tp = and kp = 

1/3 

[Z'k^y) ■ The corresponding dimensionless variables 
will be denoted with a tilde. In this way, the dependence 
on density N jV is implicit and we are left with a depen- 
dence on the inverse of the scattering length rj — 1/a and 
on the range of the potential a only. 

To begin our discussion, we show in Figjl] (a) the nor- 
malised gap A(e, (7)/Amax(o'), as a function of e = tkl 
and a = a/kp, where A,„ax = max[A(e, ct)]. Although 
these calculations were performed for the strongly inter- 
acting regime (fy = 0), the results are completely general 
in comparison with other values of the scattering length. 
We find that as long as tr < 1, the gap is a constant 
A « Amax for etr^ < 1, then abruptly drops to zero. As 
CT — ^ 0, we recover a constant gap for all e. For moderate 
and large values of ct, the gap shows a peak. As we see 
below, this novel dependence of the gap for large values 
of k dramatically changes the behavior of the thermo- 
dynamic properties. To complement the facts provided 
by the dependence of the ground state properties on the 
non-zero range potential, in Figjl] (b) we plot the behav- 
ior of the quasiparticle excitation energy Ek for several 
values of a. These results are in good agreement with 
those obtained for finite temperatures [42] . 

We now address the main point of this letter. This is 
the behavior of the contact variable C, eq. for the 
whole crossover at T = 0. In addition to its physical rele- 
vance, we discuss the contact variable because it has been 
noted recently [3D] that C is very sensitive to the contact 
approximation in the BCS mean-field approach. Fig. [2] 
(a) shows C = Ckp/ {Nep) as a function of 77 = l/{kpa) 
for different values of the dimensionless range a. First, 
we highlight the dashed black line that corresponds to 
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FIG. 2. (color online) (a) Dimensionless contact variable 
C — Ckp / (Nep) as a function of inverse scattering length 
fj — iKkpa). The color curves are the results of the full 
BCS equations for finite values of a, coded in the colored bar 
above. The dotted line is the BCS solution in the contact ap- 
proximation neglecting the HF terms. The continuous black 
lines are HYL limit, (b) The value of Tjmax where the contact 
takes its maximum value. We note that it obeys a very simple 
scaling fj 

max ~ 

l/(4a). 

the solution of the BCS equations in the contact approx- 
imation but completely neglecting the Hartree-Fock (HF) 
terms, see eq. This approximation has been deemed 
as inappropriate (30] to describe the system since in the 
deep BCS regime, 77 — —00, it predicts an exponential 
decay, as can be seen in Fig. [2] (a) [dashes] . 

The correct behavior of C is believed to be given by 
the weakly-interacting approximations of a Fermi liquid 
in the BCS side and the diatomic version of a Bose gas 
with repulsive interactions in the BEC side. Both are the 
celebrated expressions of Huang- Yang- Lee (HYL) [43]. In 
the very deep BCS and BEC regimes, those are simply 
HF corrections. In the BCS side one finds that the lead- 
ing term is C « 2/{?,'Kjf) while in the BEC side C « 2??. 
These are shown with solid black lines in Fig. [2] (a). As 
expected, the crossover cannot be described within such 
an approximation. In addition to these lines, with black 
dots, we also plot the mean-field BCS theory, in the con- 
tact approximation, with the HF terms. We see that as 
f) — >■ ±cxj, this solution approaches HYL result but, in the 
strongly interacting limit diverges. The latter behavior 
does appear to invalidate the use of mean-field theory. 
However, the divergence is due to the contact approxi- 
mation and not to the mean-field approach itself, as we 
discuss now. In color lines we show the finite range cal- 
culations. In general, we find that as a is very small, the 
contact C vanishes algebraically as f?"^, in partial agree- 
ment with HYL. The actual HYL behavior is followed 
for (7^1, while for a <^ 1 this behavior appears for 
much larger -fj. We deduce from the finite range calcu- 
lations that in the limit tr — >■ the HYL contributions 
become less important and the BCS theory with contact 
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FIG. 3. (color online) Thermodynamics at unitarity 77 = 0. 
(a) The chemical potential fi, the maximum value of the gap 
Amax, the pressure p and the contact variable C, as functions 
of the interaction range a. (b) Bertsch parameter /3 as a func- 
tion of cr. The numbers on the right correspond to references 
of different calculations. Dashed lines correspond to the pre- 
dicted values at unitarity of the BCS theory in the contact 
approximation without the Hartree-Fock terms. All quantities 
are dimensionless. 



approximation, but without HF terms is the correct re- 
sult. In the BEC side a novel behavior emerges due to 
the finite character of tr. The contact variable C peaks at 
a value that depends on a and then decays to zero again 
as 77 increases. The value of T^max at the peak scales sim- 
ply as (T~^, see Fig. [2] (b). We further note that as a 
keeps increasing, the location of the peak not only tends 
to 77 = 0, i.e. to unitarity, but the value of C becomes 
sharply peaked, suggesting that the crossover at unitar- 
ity becomes a true quantum phase transition. Although 
the transition may appear unrealistic since one needs an 
infinite interatomic range, we would like to recall that 
the present theory is a mean-field one. Such a transi- 
tion is reminiscent of the Kac limit [H] of the van der 
Waals fluid. The latter is the epitome of the equation of 
state of a real fluid, showing both first and second order 
phase transitions, yet it does not represent any real fluid. 
It remains to experimentally verify whether the present 
Kac limit is a precursor of a true phase transition for real 
atoms with finite range interactions. 

In the unitary limit (|a| — > 00) the interactions are 
strong and in the zero range limit the system posses uni- 
versal properties |45j . However, in the contact approxi- 
mation the HF terms in the Hamiltonian [last line (1)] 
become divergent at unitarity. We find that by using 
a finite range interatomic potential the terms are finite 
and are clearly not proportional to the scattering length. 
Moreover, the contact approximation can only be used 
in the weakly interacting regime. Since in our treatment, 
we never use the contact approximation, we never face 
any divergence. Thus, we find that the crossover in gen- 
eral is smooth as one would expect, except when ct ^ 
where the system approaches the Kac limit [44]. First, 
we note that, as a < 1, the exact calculation approaches 
that of the usual BCS mean-field calculation with the 
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contact approximation, but without including the diver- 
gent HF terms, dotted lines in Fig. [3] . We have that 
for small interatomic range, all thermodynamic quanti- 
ties approach the MF universality class For larger 
values of a evident deviations appear. In Fig. (b) we 
show Bertsch parameter /? = Eint/Eo — 1 [46] (Sint and 
Eq are the energies of the interacting system and non- 
interacting system respectively) for the present calcula- 
tion together with a collection of corresponding values 



for different type of models and calculations [81-fT^ [T51 - 
Ha [m EU Ea |46]. we clearly see that all those values 
are within ct < 1 of the exact calculation using the BCS 
ansatz. As one would expect, as the interaction range in- 
creases at unitarity the system starts to deviate strongly 
from the ideal Fermi gas (/3 — 0) and the result of the 
contact approximation {(3 = —0.41). Full elucidation of 
this important point is beyond the scope of this article. 

This work was partially supported by grant IN108812- 
2 DGAPA (UN AM). 
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